Extracting the Landscape and Morphology of Aging Glassy Systems 
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We propose a way to analyze the landscape geometry explored by a glassy system after a quench 
solely based on time series of energy values recorded during a simulation. Entry and exit times for 
landscape 'valleys' are defined operationally by the occurence of anomalous energy changes revealed 
by barrier and energy records. Linking these non-equilibrium events — or 'earthquakes' — to the 
record statistics of the thermal noise immediately leads to the In t relaxation behavior ubiquitous in 
glassy dynamics. Aging of Ising spin glasses in two and three dimensions is studied as a check for a 
number of low temperatures and lattice sizes. A simple picture emerges on e.g. 1) the scaling with 
temperature and system size of the energy barriers as a function of Int, 2) the scaling with system 
size of the lowest energy seen in each valley. 

PACS numbers: 02.70.Uu ; 05.40.-a ; 75.50.Lk 



Macroscopic properties of thermalizing glassy systems 
depend on the time (age) elapsed after the quench into 
the glass phase: on time scales shorter than the age, a 
state of pseudo thermal equilibrium holds locally within 
metastable regions of the landscape, while on longer ob- 
servational time scales the relaxation is manifestly non- 
stationary Some models of aging build on heuris- 
tic assumptions on the morphology of low energy exci- 
tations in real space 0, 0], while others start from a 
coarse grained state space 0, 0, IS El or > more simply, 
assume a given distribution of exit times from traps in 
configuration space [lfj . Still, two related and fundamen- 
tal questions remain unanswered: 1) which dynamical 
events mark the transition between equilibrium and non- 
equilibrium dynamics, and 2) how can the corresponding, 
dynamically selected, structures in configuration and/or 
real space be identified? 

In models of driven dissipative systems with multiple 
attractors 0, 0, El , marginally metastable attractors 
with an a priori negligible statistical weight are nev- 
ertheless typically selected by the dynamics, an effect 
which underlies memory and rejuvenation effects anal- 
ogous to those observed in the thermalization of e.g. 
spin glasses after a quench. One can speculate that 
similar mechanisms could generally be present in glassy 
systems with an extensive number of metastable attrac- 
tors. To guarantee the dynamical relevance of the at- 
tractors found and/or their real space counterparts, the 
algorithms utilized should match the dynamics in ques- 
tion, e.g. thermalization after a deep quench. This how- 
ever is not the case for standard pro cedures as e.g. the 
Stillinger- Weber approach [ij], Ha, ll3 , genetic al- 

gorithms 20] and energy minimization of excitations of 
fixed volume ElEi. 



A complementary approach is presented below which 
solely relies on statistical information collected during 
the unperturbed relaxation following an initial quench. 
The approach identifies a temporal sequence of non- 
equilibrium events, or 'earthquakes'. These entail major 



reorganizations of the landscape (1) and subdivide the 
states visited into an ordered sequence of sets, or 'val- 
leys' wherein (2) the equilibrium-like fluctuation dynam- 
ics unfolds. Both (1) and (2) are required for a consistent 
description, and both properties have been checked ex- 
tensively for short range spin glass models in 2d and 3d. 
For these models, the non-equilibrium aging of physical 
quantities as e.g. the size of the barriers overcome de- 
pends simply on the 'valley index'. This is a proxy for 
the system age marking the position of each valley within 
the ordered sequence of valleys visited. 
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FIG. 1: A fictitious time series of collected energy values illus- 
trates our sampling procedure and our definition of a valley 
(= gray area). Conspicuous events, i.e. energy and barrier 
records, are denoted by circles and squares. 

A discussion of the free energies associated with the 
fluctuation regime inside the valleys is left to forthcoming 
publications and the variable t will denote the age of the 
system, except when otherwise indicated. 

Method: A multitude of inequivalent metastable 
states produces the time inhomogeneous dynamics gen- 
erally observed in aging systems. In most cases, these 
states are sets of configurations surrounding local energy 
minima, i.e., qualitatively speaking, valleys in the en- 
ergy landscape. Since high-lying and shallow local en- 
ergy minima greatly outnumber deep minima, the initial 
quench typically produces a poor local minimum and an 
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aging system subsequently encounters a series of progres- 
sively higher energy barriers [23I I24I Ejf which delimit 
progressively deeper valleys. The values of the highest 
barrier B and lowest energy E — henceforth best so far 
or simply bsf barriers and energies — stand out as crucial 
events, and will serve to tag the valleys as they appear 
in the energy landscape. 

The simulations start from a random, T = 00, con- 
figuration and run at temperature T. Barrier values are 
computed as energy differences from the current bsf en- 
ergy state. The bsf energy and barrier are initially 00 
and zero, respectively. Each run produces a time ordered 
sequence of E's and B's as e.g. EBBEEBBBE in Fig.ffl 
Contiguous bsf energies pertain to a trajectory 'tumbling' 
down in a valley. Being mainly interested in the configu- 
ration where the tumbling stops, we only keep the lowest 
E within each group. In Fig.^the datum removed is rep- 
resented by the '(g)' symbol. Similarly, we want to avoid 
registering the numerous bsf barriers of roughly the same 
height which appear in close succession while the system 
goes through a high 'saddle'. A barrier record is therefore 
kept only if the energy experiences a sufficiently deep lo- 
cal minimum before the next global barrier maximum is 
encountered. The precise definition of 'sufficiently deep' 
involves a positive number x < 1. For x = 1, the lo- 
cal minimum must be as low as the current bsf energy, 
while for x = all bsf barriers are kept. Interestingly, 
our results are insensitive to the value of x in the range 
0.25 — 0.75. In terms of Fig. ^ we oru y keep the full 
squares and circles, thereby reducing our string even fur- 
ther to its final state: EBEBBE. 

For speed, we rely on the Waiting Time Method 
(WTM), a rejection-less Monte Carlo scheme [2^ well 
suited for problems where N variables contribute addi- 
tively to the energy through local interactions. Based on 
the local field, each variable is stochastically assigned a 
flipping time and the variable with the shortest flipping 
time is updated together with the local fields affected 
by the move. The sequence of WTM moves equals in 
probability Ejfil the sequence of accepted moves in the 
Metropolis algorithm. The current flipping time, hence- 
forth simply 'time', corresponds to a Metropolis sweep, 
and to the physical time of an experiment. 

Results: We consider L d Ising spins placed on a regu- 
lar lattice of linear size L which interact via the Edwards- 
Anderson Hamiltonian 

n{a) = -\Y, J ^ )s t ) - (D 

In this formula, s\°^ is the spin value at site i for con- 
figuration a, the couplings are symmetric in their 
two indices and vanish unless i and j belong to neighbor 
sites. In the latter case they are drawn independently 
from a Gaussian distribution of unit variance. The dy- 
namics is nearest neighbor thermal hopping with detailed 
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FIG. 2: Main figure: the average number tTb of barrier records 
versus time for three selected system sizes in 2d at T = 0.4. 
Left insert: the logarithmic slopes of Tlb as a function of the 
linear system size L for 2d (T — 0.4 and 0.6) and 3d (T = 
0.4). The upper (lower) scale refers to 2d (3d) systems. Right 
insert: the variance to average ratio for ns(t) for the data in 
the main figure. 

balance. The runs spanned over several time decades, al- 
ways bringing the system down to quite low energies, 
typically around —1.66 per spin. Large systems tend to 
have more shallow valleys than smaller system. To ob- 
tain meaningful comparisons all data originating from the 
first ten time units after the quench were discarded, and 
the value of the valley index was shifted up to one unit 
in our scaling plots in Figs. and 0] 

A basic quantity is the number of events (barrier or 
energy records) which on average occur in [0,t]. Since 
achieving a record becomes increasingly difficult, the rate 
of events will decrease in time and the dynamics decel- 
erates. Interestingly, a time homogeneous dynamics is 
restored by using the natural logarithm In t as time vari- 
able. This is shown in Fig. |2J where the average number 
ns(i) of barrier crossing events is plotted as a function 
of Int. The data shown in the main figure belong to 2d 
systems, while similar data in both 2d and 3d were omit- 
ted for clarity. The left insert shows the logarithmic rate 
of events dns(t) / dint, in 2d as well as 3d, as a function 
of L. Finally, the right insert shows the ratio of the vari- 
ance o-^ B (t) to nB(t). Both quantities were estimated 
for 3000 independent runs: The variance and the aver- 
age follow the same logarithmic fashion, with the former 
trailing the latter due to an initial lag. A similar statis- 
tics is also approached by the number of bsf records — 
and hence the number of valleys — discovered on average 
in [0,i], for sufficiently large systems. A logarithmic de- 
pendence of macroscopic averages is ubiquitous in com- 
plex metastable systems. Beside the present examples it 
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FIG. 3: Average energy barriers records B separating contigu- 
ous valleys in 3d spin glass landscapes are scaled with L and 
T, as indicated in the ordinate label and plotted versus n. All 
combinations of L = 8, 10, 12, 14, 16 and T = 0.3, 0.4, 0.5, 0.6 
and 0.7 are included, plus L — 18 and T = 0.4. 



is found in e.g. driven dissipative models fl3j ]. vibrated 
granular syste ms |281 123 ] , Lennard- Jones glass es l30| , ni- 
tridic ceramics |3l| and population dynamics |32| . 

The approximate equality between average and vari- 
ance points to the relevance of log-Poisson statistics de- 
scribing the number of magnitude records achieved in 
a sequence of t independent random numbers |l2| . The 
name stresses that lni, rather than t, is the relevant time 
variable. For data stemming from a nearest neighbor 
random walk in a landscape, log-Poisson statistics indi- 
cates a de- correlation of the energy function between con- 
secutive barrier records, and a fortiori during the time 
spent within each valley. Since the correlation time and 
the relaxation time (for the internal thermalization) can 
be generically identified |3^|, we conclude that thermal 
quasi-equilibrium is achieved within a valley. As a conse- 
quence, the morphology of quasi-equilibrium fluctuations 
around the bsf state of each valley can be isolated from 
the non-equilibrium dynamics and analyzed separately. 

The scaling of the logarithmic rate of events with L is 
shown in the left insert of Fig.[21to increase from an initial 
value of 1 to approximately 2 for large system sizes. With 
Poisson statistics, this rate is proportional to the num- 
ber of lattice 'regions' evolving independently. As the 
spins are quiescent most of the time, a record in one re- 
gion leads to an overall record as well. Thus, two length 
scales [2^ might in principle coexist: the linear size of 
the quasi-equilibrium fluctuations and the extent of the 
non-equilibrium 'earthquakes' leading from one valley to 
the next. Figs. [21 and 0] imply that the latter scale is 
asymptotically oc L, i.e. earthquakes are extensive or 
near extensive events. This is further confirmed by the 
distribution (not shown) of the Hamming distance (i.e. 
the number of spins whose orientations differ), between 
the bsf configurations of valleys n and n + 1 . The distri- 
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FIG. 4: The average difference A between the lowest ener- 
gies of the first and the n'th valley is plotted versus n. The 
ordinate is scaled with L as indicated. The system and the 
combinations of L and T shown are as in Fig. 3, except for 
an additional T = 0.2 data set (upper triangle), which falls 
outside the main trend. The remaining spread in the data col- 
lapse cannot be removed by introducing an additional scaling 
parameter, and is likely indicative of an underlying distribu- 
tion of minima. 



bution is exponential, and its average grows linearly in n 
and close to linearly in L d . 

The scaling with L and T of the energy barrier B sep- 
arating two contiguous valleys in 3d systems is shown 
in Fig. [3J All points are quenched averages over 3000 
independent runs for the L and T values specified in 
the caption. As mentioned, the data have been shifted 
by up to one unit along the abscissa, which mirrors 
the arbitrariness of the initial bsf states — conventionally 
recorded after ten time units. As indicated in the plot, 
B cx L 115 T 17 . We have not attempted a quantitative de- 
termination of the uncertainty on the scaling exponents, 
but changing the last significant digit visibly affects the 
quality of the plot. 

The above considerations apply equally well to Fig. 
which, for the same set of runs, depicts the average dif- 
ference A(l, n) between the bsf energy of valley 1 and n. 
A(l, n) is independent of T, as one should expect, and in- 
creases almost linearly with the valley index. Refs. 0>0 
analyze a mesoscopic model with this exact property. 
The curvature for large n might partly be due to a sys- 
tematic error, since not every run enters the last few 
valleys within the maximum time span allotted to the 
calculations. For very long time scales, the curve will ap- 
proach the constant value of the difference between the 
first bsf energy and the energy of the ground state. 

The strong L 2 ' 5 size dependence of A(l,n) is no- 
ticeable: since moving from one valley to the other in- 
volves the rearrangements of 0(L 3 ) spins, it seems likely 
that the set difference between two bsf energy config- 
urations be a sponge-like object of the sort discussed 
in 0, ■ The energy difference between different bsf 



4 



states is clearly not 0(1) as expected in mean-field. It 
seems unlikely that going to lower energies could radi- 
cally change the scaling of A with system size. We stress 
again that the WTM and the Metropolis algorithm being 
equivalent, everything could have been done in the con- 
ventional way, albeit with much greater numerical effort. 

In summary, the distinction between pseudo- 
equilibrium fluctuations and strongly non-equilibrium 
events in aging dynamics is emphasized. We link the 
latter events to barrier and energy records marking 
the appearance of new valleys as the landscape unfolds 
during aging. The homogeneity of the event statistics in 
In t is explained theoretically in terms of the statistics of 
records. 

The description of spin glass non-equilibrium dynam- 
ics as a log-Poisson process indicates that thermal aging 
can be treated on the same footing as driven dissipative 
dynamics 0, 0] and evolution models 0, [35| , in spite 
of obvious differences at the microscopic level. It also im- 
plies that, on each time scale, glassy dynamics typically 
selects the least stable among the metastable configura- 
tions available. While this feature has been noticed in 
more restricted contexts 0, 0, ^| , the present study 
suggests that it might provide a powerful and general 
paradigm for glassy relaxation. 
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